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ABSTRACT 

An analysis of jet stream structure and energetics 1s performed 
utilizing an objective scheme that calculates wind camponents nomnal 
and parallel to the jet axis, with provisions to explicitly calculate 
the eddy components of variables. Most of the meteorological data used 
is generated by a diagnostic-balanced, w-equation model. Results indi- 
cate that above 500 mb observed winds and temperatures should be used 
versus computer generated data. 

There is same imication that between the jet entrance and exit 
regions the mean transverse vertical circulation reverses from thermally 
direct to thermally indirect. Kinetic energy is imported from below 
the jet into the jet stream region primarily by the mean pressure 
interaction at 500 mb. 

Thermal balance calculations indicate that the mean horizontal 
temperature advection and vertical transverse motion are primarily 


responsible for maintenance of the temperature structure about a jet. 
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I. INTRODUCTION 

Since the mid-1940's, when the existence of jet streams in the 
troposphere was documented by Rossby and Staff (1947), many studies of 
the Polar Front Jet Stream have been conducted. Notably, these studies 
have yielded a considerable amount of empirical knowledge about the 
Polar Front Jet. Unfortunately, the energetics of this jet are less 
well understood. 

One important reason for this lack of understanding has been 
pointed out by Mahiman (1970). Considerable difficulty exists 7 
defining an adequate curvilinear coordinate system for such a strongly 
time-dependent phenomenon. Some investigators have attempted to over-: 
come this difficulty and have conducted studies of jet streams in 
curvilinear coordinates, with most interesting results. Riehl and 
Fultz (1957, 1958) in their steady, three-wave rotating dishpan experi- 
ment, investigated the dishpan jet stream in curvilinear coordinates, 
centered on the jet axis. They found net ascent poleward and net 
descent equatorward of the strong westerlies when averaging was per- 
formed around latitude circles. However, when the averaging was per- 
formed along the jet axis, the vertical eto was reversed. 

This implied that the mean transverse circulation about the dishpan 
jet was thermally direct, while the mean meridional circulation in 
the same latitude operated in a thermally indirect sense. This 
suggested that the mechanisms responsible for maintaining the jet 
streams in a geographic coordinate system may be different from those 


in a curvilinear coordinate system. 





The most detailed study of an atmospheric jet stream in curvilinear 
coordinates was performed by Krishnamurti (1961), on the subtropical 
jet stream, covering the period December 1955 and January-February 
1956. The subtropical jet stream is particularly suitable to analysis 
in curvilinear coordinates since its day-to-day variability is rather 
small, especially when compared with the polar front jet. Krishna- 
murti's study also showed the existence of a thermally direct trans- 
verse circulation about the jet core. This agreed with the net 
transverse circulation about the jet in curvilinear coordinates.. The 
mean meridional circulation in the latitudes of the subtropical jet 
operates in the same sense as the transverse circulation, in contrast 
to the dishpan jet. 

The preceding studies avoided the time-dependent problem associated 
with jet streams. In the case of the dishpan jet the flow was in 
steady state with only longitudinal translation, and in the subtropical 
jet case the flow patterns were nearly stationary. In addition, both 
studies were limited to same extent, in that no explicit calculations 
of eddies in a curvilinear coordinate system were made. 

Until recently no attempt has been made to extensively investigate 
the Polar Front Jet Stream in curvilinear coordinates since it is the 
most time-dependent of the major tropospheric jets. Mahiman (1970) 
minimized this time-dependent problem by choosing a Polar Front Jet 
where the time variations were not particularly large. While his 
study was quite extensive, it was limited primarily by two factors. 
First, an excessive amount of time was required to define the grid 
system, extract data, and perform the necessary calculations in the 


curvilinear coordinates. The second difficulty was in explicitly 
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computing eddy components of the wind in a curvilinear coordinate systen. 
In this system an eddy quantity was defined as the value of a variable 
at a given point minus its mean value along the same line parallel to 
the jet axis. As Mahiman (1970) pointed out, in a curvilinear coordi- 
nate system, the system tends to be aligned nearly parallel to the flo. 
Thus, the computation of wind eddies was extremely sensitive to the 
orientation of the jet axis and small variations in the wind flow. 

The first intent of this study was to develop a method for rapid, 
objective analysis of a jet stream. Meteorological input data was 
generated by a diagnostic-balanced, w-equation program (Krishnamurti 
1968), which required only geopotential heights. Thus both the 
divergent and nondivergent components of the wind are included in the 
study. Some discussion of the advantages and disadvantages of the data 
used, plus some steps necessary to minimize the disadvantages, is 
contained in Section II-B. 

The objective computational scheme defined a coordinate along a 
jet stream axis and coordinate lines at distances 2.5 and 5.0 degrees 
latitude normal to the jet axis on either side (Fig. 1). Then the 
shceme linearly interpolated data from the various fields; w,T,%,U and 
V camponents of the total balanced winds. The scheme calculated wind 
components normal and parallel to the jet axis. Line or area averaging, 
where appropriate, was performed to calculate the eddy components of 
variables. A detailed description of this scheme is contained in 
Section II-D. 

Once the above computational scheme was developed, several checks 
were performed. First, results of the data interpolated fram the 


initial data fields by the camputational scheme were compared with 
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independent interpolations done by hand. Second, using the same time 
periods and jet axes at 300 mb that Mahiman (1970) used, the jet stream 
kinetic energy balance was calculated. The basic data fields and cam 
putational procedures were discussed by Mahlman (1970). The results 
obtained were then compared with Mahlman's results based on subjective 
interpolations and hand calculations. 

The results of these checks were quite successful. Comparison of 
the data interpolated from the initial data fields by the computer and 
by hand showed very small differences, especially when averaged over 
the five time periods involved. This meant that the data interpolated 
by the objective scheme could be used with confidence in calculating 
the jet stream kinetic energy balance. Camparison of Mahlman's and the 
objective scheme calculations for the jet stream kinetic energy 
balance also showed good results which will be presented in Section IV. 

Having satisfactorily developed the objective computational scheme, 
three further objectives for this research were set. First, due to 
certain inadequecies in the available data, which will be discussed in 
detail in Section II-B, jet axes at 200 mb were defined. The jet 
stream kinetic energy balance with respect to these axes was then cal- 
culated as it was for the axes at 300 mb mentioned previously. Second, 
Since Mahlman's results were based on five time-period averages, it 
was felt that it would be of some interest to study a single time 
period. Therefore, the jet stream for 16 November 1200 GT, 1966, 
was analyzed in detail. This period appeared to present the most 
explicit situation of a jet stream maximum of the five time periods 
studied. 
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Lastly, Mahlman (1970) discussed the role of eddy heat flux in the 
maintenance of the temperature structure about the jet. Mahlman's 
analysis indicated that eddy heat flux might counteract the expansion 
cooling and compression heating caused by the mean circulation about 
the jet. Therefore, the third objective of this research was to 
explicitly investigate the role of eddy heat flux in the maintenance 


of the temperature structure about a jet stream. 


II. RESEARCH PROCEDURES 

A. DESCRIPTION OF THE JET STREAM COORDINATE SYSTEM 

The jet stream coordinate system for this project was defined, 
consistent with Mahlman (1970), in the following manner. A central 
coordinate was established, in this case, along the jet stream axis. 
Parallel coordinate lines were defined at distances 2.5 amd 5.0 degrees 
latitude, normal to the jet axis in either direction. In an (s, n, p) 
curvilinear coordinate system, s and n were respectively the coordinates 
parallel and normal to the jet axis, and p was taken as the vertical 
coordinate. Separation between grid points along the s coordinates, 
other than the jet axis, was not uniform. This was the major difference 
between this grid and the one Mahlman used with unifomm separation along 
all s coordinates. Once this jet stream-oriented coordinate system was 
defined by specifying an axis at a single pressure level, it is consistent 
through all levels in the vertical (See Fig. 1). 
B. METEOROLOGICAL DATA USED 

The case selected for analysis was an area centered over continental 
United States as described by Mahlman (1970). Data at five synoptic 


observation times was used, from 15 Novenber 1200 GMT to 17 November 
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Curvilinear grid for jet stream of 15 November 1966,. 
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Fig. 2 Geopotential heichts(m) at 300 mb on 15 November 
1966, 1200 GNT 
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1200 Gyr, 1966. During this period the flow at jet stream level was 
anticyclonic between troughs off the west coast and over Hudson Bay in 
the east. Although there was some change in the location of the jet 
stream, no significant development or large-scale system mene 
occurred. The jet stream axes did show a tendency to shift to the 
southeast with the jet maxima traveling eastward along the axis. Thus, 
this case appeared to present a satisfactory steadiness of the jet 
axis for analysis (See Fig. 2). 

As mentioned in the introduction, all meteorological data in this 
research were generated by the diagnostic balance, w-equation model, 
developed originally by Krishnamurti (1966, 1968). In this model, the 
equations are scaled for motions of characteristic Rossby number less 
than unit (Ro<i). Initial input data was obtained from the National 
Meteorological Center's objectively analyzed geopotential fields at 
1000, 850, 700, 500, 300, and 200 mb. The stream function () was 


related to the geopotential field through the non-linear balance 


equation, 

Ve£VW+27 Cees =V26 (1) 
If the ellipticity condition for Eqn. (1) 

V26+1/2£2-VE*Vy>o (2) 
is not satisfied, then Egn. (1) is replaced by 


Ve£VYHV 26-27 (u,1V5) (3) 
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The meteorological data output was in staggered levels as shown in 
Fig. (5). Fields of , %, and <~ appear at the 1000, 800, 600, 400, 
and 200 mb levels, while w, T, and © are given at the 900, 700, 500, 
300 mb levels. The horizontal grid points were spaced 2.5 degrees 
jatitude and longitude apart. The grid extends from 70W to 135W and 
from 25N to 60N. Six additional grid points in the zonal direction did 
not contain initial data, but were used to provide cyclic continuity 
for any given dependent variable. These grid point values were inter- 
polated fram values at the first two and last two of the grid points 
with real initial data. There was no staggering of variables in the 
horizontal. For a concise description of the computational procedures 
of this model, see Krishnamurti (1968) and Mahlman (1970). 

Some limitations of the data have appeared. First, the jet stream 
chosen, while of moderate intensity, was not particularly well-defined. 
For most of the time periods, two distinct jet axes, whose separation 
was less than 10 degrees latitude, were present. During the last two 
time periods, 17 November 0000 GMT and 17 November 1200 Gv, 1966, a 
moderate jet stream associated with the low center over Hudson Bay 
impinged on the northeast corner of the grid, merging with the jet 
stream being analyzed. This caused same difficulty in defining axes 
and in interpreting the results. The manner in which this problem was 
handled is discussed in Section II-<. 

Second, the ellipticity problem of the balance equation along with 
erroneous geopotential heights adjustments had to be considered. When 
the ellipticity condition for Egn. (1) was not satisfied, Eq. (3) was 


substituted for Eqn. (1). This adjustment was frequently required on 
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the anticyclonic side of the jet stream. Ellsaesser (1968) has com- 
pared the calculated winds of the balance equation with the observed 
winds. In places he found the mean jet stream at 300 mb was seventy 
percent stronger than indicated by the balance equation. Calculated 
Wind speeds just equatorward of the jet stream were also found to be 
Significantly lower. Calculations bear out this disadvantage of the 
non-linearized balance equation. 

Errors in data handling also appeared. This problem was not 
particularly serious, since it was obvious when the error was large. 
Subjective adjustments to the geopotential heights were made when 
necessary. 

Since it was not an objective of this research to analyze the 
balance equation model, no attempt was made to change the model to 
minimize the ellipticity problem. However, when interpreting the 
results this limitation was kept in mind. 

Third, having the vertical motion defined at different levels than 
the kinetic energy and temperature is an obvious limitation. The 
vertical advection terms were calculated by linearly interpolating the 
values between levels. While not considered serious this problem must 
be kept in mind when interpreting the results. 

C. DESCRIPTION OF JET STREAM AXES ANALYZED 

For descriptive purposes the jet axes analyzed were divided into 
two groups, those based on observed winds and those based on the com- 
puted (balanced) winds. The first group used jet axes defined fran the 
National Meteorological Center's maps of observed 300 mb winds and 


those were the ones used to campare with Mahlman's results. 
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Fig. 3 Jet stream axes at 300 mb, defined from 
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Fig, 4 Jet stream axes at 200 mb, defined from 
the total balanced winds 
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The second group used jet axes at 200 mb defined fram balanced | 
winds (See Fig. 4). The jet streams thus specified were of moderate 
strength, fairly well-defined, and presented a reasonable representation 
of the observed winds at 200 mb. Unfortunately, the definition of the 
jet stream axes was not particularly good east of 95W. Consequently, 
only that portion of the jet axes west of 95W was analyzed. These five 
axes gave a good picture of jet stream flow fram SW to NE but only a 
half wave-length was being analyzed. 

D. DESCRIPTION OF OBJECTIVE COMPUTATIONAL SCHEME 

A computer program was developed for the IEM 360/67 system of the 
W. R. Church Computer Center which performed the calculations described 
in the following section. The remainder of this section is a general 
description of the canputational scheme used. Primary emphasis is 
placed on operations performed and principles applied. The generality 
and flexibility of this program is also discussed in some detail. For 
a detailed discussion of the mechanics concerning how the operations 
were performed, see Appendix A. 

The computer computational scheme was divided into four subroutines: 

1. data input 

2. development of the curvilinear grid to be superimposed on 

the geographic data fields at all levels 

3. interpolation of data from geographic to curvilinear system 

4. computations performed using a curvilinear coordinate system. 
The data input subroutine was simply the reading of the initial data 
fields and the definition of the jet stream axis in geographic 


coordinates. 


1b) 








The curvilinear coordinate system developed in the next subroutine 
was a simple design (See Fig. 1). The points along the jet stream axis 
were located at two degree latitude intervals along the jet axis. Then, 
utilizing the slope of the jet axis at each point in a Scher described 
in Appendix A, grid points were established at 2.5 and 5.0 degrees 
latitude normal to the axis. The grid was constant for all levels in 
the vertical. The area enclosed by the curvilinear grid and the distance 
between grid points was also calculated. 

Data was interpolated fram the geographic grid to the curvilinear 
grid in the third subroutine. A linear interpolation in two-dimensions 
was used for this transformation. After the data was interpolated fran 
geographic to curvilinear coordinate points, a linear interpolation in 


one-dimension (vertical) was performed to obtain values of w and 6 


between the staggered levels. 

In the fourth subroutine with fields of temperature, vertical 
motion, geopotential heights, and the u and v components of the total 
balanced wind at the curvilinear grid points available, the first cal- 
culation was to convert the u- and v wind components to components 
parallel and normal to the jet axis at each point. This is done by the 
scheme described in Appendix A, utilizing the jet axis slope. Next, the 
eddy components of each variable were calculated. Either a line or area 
averaging scheme, as appropriate, was used to obtain the mean of a 
variable along the coordinate. Then this mean was subtracted from the 
value of the variable at each point. Both eddy and mean values were 
required for the kinetic energy balance and the eddy heat flux calcu- 
lations. To facilitate the discussion, the equations used and the 


results are discussed in Sections III and IV. 
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One of the advantages of this program was its generality and flexi- 
bility. This program was made general so that it might be used with 
other than polar-stereographic maps and 27 by 15 grid arrays. The pro- 
gram was flexible in that there is no restriction that the 2 coordinates 
be at 2.5 and 5.0 degrees latitude from that jet axis, or that there be 
only two coordinate lines on either side of the jet axis. Thes 
coordinate lines were not limited to 2 degrees latitude distance along 
the axis. There may be any number of s coordinates, and of any interval 


length desired. 


IIIT. EQUATIONS 


In this report, the averaging notation was that used by Mahliman (1970). 


The mean of an arbitrary quantity, Q, was given by [Q] (s)° The sub- 
script enclosed by small parentheses indicates the coordinates over 
which the averaging was performed. The deviation of Q from that average 
was given by (Q), i.e., Q[Q] 8) + (Q) (s) for averaging processes along 
the jet. When averaging normal to the jet, n was substituted for s. 
Consistent with Mahlman (1970) the following equation was used to 


obtain the kinetic energy balance of the jet stream. 
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The separate terms, labeled (a) through (j) have the following 
interpretation: 


(a) flux of "mean" kinetic energy through the side boundaries; 

(b), (c) flux of "mean" kinetic energy through the lower and 
upper boundaries, respectively; 

(a) conversion of "eddy" kinetic energy to "mean" kinetic energy; 

(e) an additional term in the conversion fran "eddy" to "mean" 
kinetic energy which arises because the Coriolis parameter (f) 
varies along the axis over which the initial average is 
performed. It is interpreted as a conversion term because it 
appears with an opposite algebraic sign in the "eddy" kinetic 
energy equation; 

(£) mean pressure interaction term at the side boundaries arising 
from energy conversion in an open system; 

(ig), (h) mean pressure interaction terms at the lower and upper 


boundaries, respectively; 
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(1) conversion of potential energy into "mean" kinetic energy due 
to mean transverse circulation about the jet core; 
(j) dissipation of "mean" kinetic energy due to friction and sub- 
- grid-scale mixing. 
Mahiman specified that the "mean" kinetic energy balance in the symbolism 
implied in Egn. (4) was: 


Staa,+b,+¢,+f+gth+i - "dissipation: 


"Dissipation" was defined to be a cambination of residual tems 
aor bo, Coy and j. The subscripts 1 and 2 referred to the mean and 
eddy values of a term respectively. 

To investigate the role of eddy heat flux, the first law of therm- 
dynamics was averaged along the jet axis and separated into mean and 


eddy components. 
| (1) (2) 


ne) (s) otf) (s) 9 


—se = > N53) 5) Tae 7 i M5) Gey Mg) os) 
(3) (4) 
; Pa oe) s, oo a : Sp (s) 
(5) (6) 
a o{[T 
5 ap so) a ia gh oe 


When discussing the results of calculations using the above 
equations, the alphabetic or numeric designations of terms are used. 
This is to simplify tables presented and to not have to continually 


repeat the meaning of temns. 
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IV. RESULTS 

To facilitate discussion, the results of this research are presented 
in the following manner. First, section IV-A presents the results of 
the comparision between data interpolated from the initial data field 
by the computational scheme versus independent interpolation by hand. 
Second, sections IV-B and IV-C present the results of the calculation 
of the jet stream kinetic energy balance for the previously described 
axes at 300 mb and 200 mb (see Section II-C). Section IV-B also contains 
Mahlman's (1970) energy calculation for the 300 mb jet axes 
with the results of this research. Third, in section IV-D the analysis 
of the jet stream defined at 200 mb by the total balanced winds for 16 
November 1200 GMT is presented in same detail. EBoth jet stream struc- 
ture and calculations of kinetic energy are presented. Lastly, in 
section IV-E results are presented of the investigation to explicitly 
determine the role of eddy heat flux in the maintenance of the temperature 
structure about a jet stream. The conclusions drawn fram these results 
are presented in section V. 
A. COMPARISON OF INTERPOLATION METHODS 

The results on the interpolation of the initial data fields fram 
geographic to curvilinear grid points by both the objective and manual 
techniques are shown in Figs. 5, 6, and 7. Fig. 5 shows the mean 
vertical motion relative to the jet axis averaged over the five time 
periods. The number above each point was obtained by the objective 
scheme, the value in parentheses was obtained by hand interpolation. 
No difference was greater than 10 x 1075 mb sec™!, and at more than 


50 percent of the locations the difference was less than 5 x 107° mb 


25 











Fig. 5 Comparison of mean vertical motion (Lw]J. , 10 
mb sec™’ ) interpolated by the objective “Scheme 
| versus hand interpolation 
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Fig..6 Mean wind (knots) and temperature (deg K) 
calculated by the objective computational scheme 
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Fig. 7 Mean wind (knots) and temperature (deg K) along the 
polar front jet stream averaged from 15 November, 
1200 GMT to 17 November 1200 GMT, 1966. 
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sec”!, In the mean, the vertical motion field did not agree with Riehl 
and Fultz's (1957, 1958) finding of ascent equatorward and descent pole- 
ward of the jet axis (Section I). The first and last time periods did 
show ascent equatorward and descent poleward, but these results were 
cancelled out by ascending motion poleward of the jet axis during the 
other three time periods. 

This difference in mean values is explained by the orientation of 
the jet axes used. Fig. 3 shows that the axes used did not cover an 
entire wavelength. They were predaminantly in the region of a trough 
and of a ridge. This corresponds to regions of southerly flow ea 
anticyclonic curvature. Riehl and Fultz (1957, 1958) found that for 
these specific regions, ascent occupied the area. Since these two 
regions were predominant in the averaging performed, ascent should be 
found throughout, being stronger equatorward. The ascending motion was 
somewhat stronger than Mahlman (1970) showed, but the patterns were 
the same. In both cases, ascending motion through the jet core was 
found. 

Fig. 6 and 7 show the results for the mean wind (U5) and temperature 
computed by the objective scheme and by Mahliman (1970) respectively. 
The mean winds were similar both in strength and pattern, and the mean 
temperatures were almost identical. Thus the results of this section 
indicate that the objective scheme gives comparable one to the manual 
scheme for interpolating data fram geographic to curvilinear coordinate 
points. 

B. KINETIC ENERGY CALCULATIONS FOR 300 MB JET AXES 

The results for the 300 mb jet stream axes kinetic energy calcu- 

lations for both the objective schere and Mahlman (1970) are shown on 


Table 1. The calculations were performed from 100-900 mb over the area 
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defined by the curvilinear grid described in Section II. To acquire 
further insight, and enable more complete comparison with Mahiman's 

results, the volume was divided into regions below (300-900 mb) and 

above. (100-300 mb) the jet. In addition, calculations of terms 


Ass b d and e of Eon. (4) are included at the bottom of the 


oD 
table. This enabled some estimate to be made of eddy terms which 
Mahiman included as "dissipation". Although the results with the 
objective scheme were generally comparable with the hand-calculated 
values, there were som significant differences. 

This was expected since the objective scheme used balanced winds 
and Mahiman used observed winds for kinetic energy computations. Where 
no horizontal wind calculations were involved, terms g and h, Egn. (4), 
the results were almost identical. However, a definite difference 
existed, especially at upper levels, between terms where horizontal 


wind calculations were made, t is, terms —, ays b anf f in 


aoa 
Eqn. (4). Thus, the problems converning balanced winds discussed in 
section II-B must be considered when comparing these terms. The local 
change term ( a ) best illustrated this problem, since it involved 
only us and v.. Table 1 shows that for all layers the objective scheme 
results for enn were lower than those calculated using observed winds. 


dt 
For this case this was in agreement with Ellsaesser's (1968) finding that 
balanced winds generally were lower than observed winds. 
The role of the flux of mean kinetic energy and pressure interaction 
were the most significant differences between the two methods of cal- 
culations. In the 100-300 mb layer Mahlman's results show a net gain 


of +5.7 ergs an2mb- ‘sec 7 versus the objective schemes net loss of 
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2) ecies cm7’mb7 Jsec™1£or the flux of kinetic energy, terms a, rb), and 
Cy in Eqn. (4). The pressure interaction terms (f, g, and h) in 
Eqn (4), both show a net gain, +6.3 versus +13.4 ergs cm7~’mb7!sec~ 1}. 

For the 300-900 mb layer, however, the net results are Wei close. 
The flux terms a, by, and Cyr Egn. (4) show a net loss in both cases, 
-2.9 versus -3.0 ergs an-*mb7/sec7!. The pressure interaction terms 
also showed almost identical net losses, -4.4 versus ~-4.2 ergs 
em72mb7 ‘sec™ }. 

The large difference in the upper layer 100-300 mb, had a biasing 
effect on the total layer (100-900 mb). For this layer the flux of 
kinetic energy was of the same sign, that is, a net loss, but the 
objective scheme's results were larger, -0.8 versus -3.0 ergs 
com72nb~ ‘sec 1. For the pressure interaction effects, the objective 
scheme shows a net gain, +0.3 ergs cn mb sec versus a net loss 
shown by Mahlman, -1.6 ergs cm”-’mb74sec™}. Since terms b), C), g and 
h, Egn. (4), gave essentially the same results, the terms a, and f, 
Egn. (4), in the 100-300 mb layer were the primary difference. 

These discrepancies were felt to be the result of the following 
differences. First, the objective scheme used balanced versus the 
observed winds used by Mahlman. Second, the objective scheme explicitly 
calculated V5: Thus the calculations performed were affected by the 
jet axes to the south which effected this region (Section II-B). This 
effect was not considered by Mahlman (1970). 

The results of the objective scheme and Mahlman's are comparable, 
if the problems mentioned above are considered. However, it does appear 
that the results of calculations using balanced winds at or near jet 


stream level are questionable. 
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The lower portion of Table 1 shows the calculation of eddy temms 
Aor bo Cor d, e and thus the implied value of dissipation. Calculation 
of these terms primarily involved using Us and v_, with the previously 
Peniioned problems involved. Taking these problems into account, two 
findings concerning net results appear to be significant. First, the 
net value of the eddy terms is large enough to imply that these terms, 
Ass bo, Cor d and e in Eq. (4), do have an effect on the kinetic 
energy balance. Second, the residual term j in Egn. (4) is of sufficient 
Size to suspect that dissipation of mean kinetic energy due to friction 
and sub-grid-scale mixing does occur. 
C. KINETIC ENERGY CALCULATIONS FOR 200 MB JET AXES 

The results of the kinetic energy calculations using jet axes at 
200 mb are shown in Table 2. As in the previous section, results are 
presented for a total layer and two sub-layers. The sub-layers, however, 
are different from those shown in Table 1. There were two basic reasons 
for this. First, the jet axes were at different levels, 200 mb versus 
300 mb. Since the balance wind data was available only to 200 mb, it 
was impossible to calculate the wind in the region above 200 mb. Second, 
from the results of the comparison in section IV-B and the problems 
involved with the data (Section II-B), it was felt a different division 
of sub-layers was called for. Thus, the total layer was divided into 
two equal sub-layers, 500-900 mb and 100-500 mb. This took into account 
the apparent fact that balanced wind calculations in the vicinity of the 
jet stream were less reliable. 

As shown in Table 2, one of the sifnificant results was the large 
value of the flux of mean kinetic energy through the side boundaries, 


tem a,, Eqn. (4), for the upper layer, 100-500 mb. It was large enough 
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TABLE 2 


All units expressed in erg cm mb sec 
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to completely overshadow the net results both for that layer and the total 
layer, 100-900 mb. In the lower layer, 500-900 mb this dominance by 

term ays Eqn. (4) did not exist. In this layer the mean flux of kinetic 
energy through the upper boundary, term Cc), Eqn. (4) and the eddy flux 

of kinetic energy through the side boundaries, term Ass Eon. (4), had 

the dominant roles. For this layer, the mean and eddy terms almost 
cancelled, leaving a net -1.9 ergs cm7~2mb7 4sec™ 1. 

Results for the 200 mb jet axes were not directly comparable with 
the 300 mb axis. The balanced wind output tended to merge the two axes 
into a single axis closer to the northern axis in the first three time 
periods, and the southern axis in the last two periods. Furthermore, 
they covered only a half-wavelength. However, these results did allow 
for a check on the results presented in section IV-D of the jet stream 
at 200 mb for 16 November 1200 Gr. 

D. ANALYSIS OF AN INDIVIDUAL JET STREAM AXIS 

The results of the individual jet stream analysis are presented for 
the entrance and exit regions and are further sub-divided into a total 
layer and two sub-layers, 500-900 mb and 100-500 mb. 

Fig. 8, 9, and 10 show the averaged vertical motion, temperature, 
and wind fields of both regions respectively. The results of the kinetic 
energy calculations for both regions are shown in Table 3. Table 3 also 
shows the result of averaging along the entire jet for use in making the 
camparison mentioned in section IV-C. 

Considering the total jet structure first, vertical motion in the 
entrance region was ascending equatorward and descending poleward of 
the jet axis (Fig. 8). The normal wim, V_, was equatorward below the 


jet defined at 200 mb (Fig. 10). Unfortunately, wind fields were not 
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for November 1200 GMT 
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available in the balance model above the jet stream level to check if 
the transverse winds there were poleward, thus completing a thermally 
direct circulation. In this region the nommal winds, V5. were 
considerably stronger than those found in the exit region. 

The vertical motion in the exit region was ascending throughout, 
although much stronger poleward of the axis (Fig. 8). The normal wind 
field, V., did show same poleward flow beneath the jet, but the values 
were quite small (Fig. 10). Thus, an indication of a thermally indirect 
circulation did exist, but was not particularly well-defined. The large 
normal wind values south of the jet axis at jet stream level (Fig. 10), 
are attributed to the location of a jet stream to the south (see Section 
II-B and Section IV-B). 

The temperature structure (Fig. 9), of the two regions is interesting. 
At jet stream levels, the entrance region is 3-5 K warmer than the exit 
region though the temperature patterns are the same. At 700 mb the 
temperatures for the two regions are almost identical. The large 
difference in temperatures poleward of the axis at 900 mb is attributed 
to a cold surface outbreak that occurred during this time. 

As shown in Table 3, one of the most significant results again is 
the extremely large value, for the 100-500 mb layer, of the mean and 
eddy flux of kinetic energy through the side boundaries, terms ay and 
ot Eon. (4), especially term Ao. Although large values had been found 
before, none had been comparable to tem Ao: Egn. (4), as shown in 
Table 3. These results were too large to be considered reliable, and 
could not be explained by the data problems discussed in Section II-B 


and IV-B 8 
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TABLE 3 


All units expressed in erg cm mb _ sec 
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The fact that only a portion of a wavelength was analyzed, and 
that it was further sub-divided into entrance and exti regions, must be 
considered. The findings indicate results for a segment of a jet azis, 
and not the averaged results along an tines oe The change of 
Signs between regions, Table 3, indicates balancing between regions 
along the jet axis. Since in this case the wind acceleration is much 
stronger in the jet entrance than the deceleration in the exit zone at 
jet stream level, the regions should not be expected to completely 
balance. Notice that at lower levels, where acceleration and 
deceleration are about equal, the two regions almost balance. 

Taking the above into consideration, two conclusions result. First, 
the lower layer, 500-900 mb, is below the area where shears are large, 
and the results of the objective scheme may be considered reliable. 
Second, in the 100-500 mb layer terms g, h, and i, Egn. (4) do not 
contain horizontal wind calculations, and term f, Eqn. (4) only uses Cy 
THus, some comparison of these terms between layers is possible. 

In Table 3 the entire jet results are slightly different fram the 
average of the two regions because of a 2° separation between entrance 
and exit regions caused when the grid was developed by the objective 
scheme. Notice that the entire jet results for the 500-900 mb layer 
are very close to those shown in Table 2, Section IV-C, for this layer. 
Since these results are for two different situations, it is reasonable 
to assume that the results for the entrance and exit regions are reliable. 

Comparing the 500-900 mb layer jet entrance versus the jet exit 
region, it is interesting to note the change of signs. For every term 


calculated in Egn. (4), except c,, the algebraic sign changed between 


1 
regions. In addition, the entrance region showed a net loss, while the exit 


region showed a net gain for all terms (-16.8 vs +15.3 ergs cm-*mb~ 'sec—}). 
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In the jet entrance region, kinetic energy was lost primarily by 
mean pressure interaction through 500 mb and the eddy flux of kinetic 
energy through the side boundaries. The exit region's kinetic energy 
gain was also primarily caused by mean pressure interaction through 
500 mb and the eddy flux of kinetic energy through the sides. But 
viewing the entire jet stream the upward transport by the mean pressure 
interaction in the entrance region is offset by downward transport in 
the exit region. 

E. INVESTIGATION OF EDDY HEAT FLUX 

This investigation was conducted on the same jet stream axis used 
in Section IV-D, and was divided into the same regions, jet entrance 
and exit zones. Fig. 11 shows the magnitudes of the terms in the heat 
budget equation, Eqn. (5), for both regions and at various levels. 
Since temperatures above 300 mb were not available from the balance 
model, only results for the region below 300 mb are shown. The vertical 
convergence of eddy heat flux and eddy adiabatic compression, term 5, 
Egn. (5) was negligible and not included in Fig. ll. The diabatic 
temperature change, tem 4, Eqn. (5) was not calculated, although same 
estimates of it were made and are discussed later in this section. 

Mahlman (1970) suggested that the horizontal convergence of eddy 
heat flux, term 2, Eqn. (5) might counteract the expansion cooling and 
compression heating caused by the mean circulation about the jet. 

Fig. 11 shows that except for the region poleward of the jet axis at 
400 and 600 mb, eddy heat flux does not have an appreciable effect. In 
the entrance region eddy heat flux was the smallest of the four terms 


and showed a slight net transfer of heat poleward. 
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In the exit region the eddy heat flux was camparable in magnitude to 
the mean horizontal transverse advection, term 1, Eqn. (5), at both 
400 and 600 mb. However, at these levels, both terms were small. 

In the jet entrance region Fig. 11 shows that the mean horizontal 
advection terms (1 and 3, Eqn. (5)) were almost balanced. The cold 
advection of the mean horizontal normal advection, tem 1, Eqn. (5), 
being slightly stronger than the warm advection of the mean horizontal 
parallel advection, texmm 3, Egn. (5). In the jet exit region, the warm 
advection of the mean horizontal advection, term 3, Egn. (5) was. 
approximately balanced by the expansion cooling of the mean transverse 
vertical motion, especially poleward of and at the jet axis. In the 
jet entrance this balancing of the mean horizontal advection leaves 
the mean transverse vertical motion as primarily responsible for the 
warming poleward of the jet axis in the upper levels. The mean trans- 
werse vertical motion causes cooling equatorward of the jet. In the 
jet exit region, all the terms approximately balance except at 600 and 
800 mb along the equatorward boundary, where cooling is shown. 

An interesting point shown in Fig. 11 is the reversal of roles 
between jet regions by the horizontal normal advection and the vertical 
transverse motion. Notice that in the jet entrance, the mean trans- 
verse vertical motion pattern is approximately the same as the mean 
horizontal parallel advection. In the jet exit region the mean 
horizontal normal advection pattern closely follows the mean horizontal 
parallel advection. 

The magnitudes of all terms shown in Fig. 11 decrease equatorward, 
except for the compression heating or expansion cooling due to mean 


transverse vertical motion, term 6, Egn. (5), in the entrance region. 
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d[T 

(5) 
ot 

is almost constant for all three levels across the jet. Thus, it is 


However, aS Fig. 12 shows, the local change of temperature, 


implied that to maintain this structure, the magnitude of the diabatic 
temperature change must increase equatorward. Fig. 13 shows that the 
estimated magnitude of the diabatic temperature change does increase 
equatorward, balancing the decrease equatorward of the other terms in 
em. (5) . 

Confidence in the actual values of these results was limited. The 
temperatures used were computer generated and then interpolated to the 
same levels as the balanced winds. However, fram the relative magnitude 
obtained, eddy heat flux did appear to play a role poleward of the jet 
axis. This agrees with Krishnamurti's (1961 b) findings for the sub- 
tropical jet that north of the jet axis there is an indication that heat 
flux is carried largely by daily eddies. The indication of this research 
is that the eddies do play a role in the heat flux, but not a dominant 
one. In addition, it appears that the diabatic temperature change 


equatorward of the jet axis plays a significant role. 
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for 800, 600, and 400 mb 


= 
Gi 
aie 2 , se = C00 
| ee pe: 
- Ee ee : 
6 ‘is <s 
if 
aie 
ere 
4 
is ee 
NM 5.0 ais ah aye See 5S) 


Fig. 13 Magnitude of the estimated diabatic temperature 
change, ean. (5) (in°Chr’ ) for 16 November 
1200 GMT, for 800, 600, and 400 mb. 
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V. CONCLUSIONS 

Five valid conclusions are to be drawn from this research. First, 
the objective scheme developed did provide for a rapid and effective 
analysis of a jet stream. Wind computations in the vicinity of large 
wind shears may, however, require a finer grid than used. It is 
recommended that if wind computations are desired in large wind shear 
areas, the grid be small enough to resolve these large shears. 

Second, use of total balanced winds in the vicinity of the jet 
stream are questionable. It is recommended that observed winds be 
used above 500 mb. 

Third, there was a strong indication that the transverse circulation 
was thermally direct in the entrance zone and thermally indirect in 
the exit zone of a jet maximm. This was in agreement with the 
conclusions of Murray and Daniels (1953). 

Fourth, kinetic energy is imported and exported from below the jet 
into the jet stream region primarily by the pressure interaction at 
500 mb. This export and import along the axis is not seen if averaging 
is performed along the entire length of the axis. It is suggested that 
when studying jet stream energetics, the jet entrance and exit zones be 
included in addition to calculations along the entire length. It 
appears that some important features in maintaining a jet maximum tend 
to be balanced between these two regions. 

Fifth, there was an indication that mean horizontal temperature 
advection and mean transverse vertical motion paly a significant role 


in the maintenance of the temperature structure poleward of the jet 
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axis. Equatorward of the jet axis there was a strong indication, fron 
estimated results, that the diabatic temperature change plays a 


Significant role in maintaining the thermal structure. 


48 





APPENDIX A 

Throughout this report frequent reference to an objective compu- 
tational scheme was made. A general description of. the is was 
given in Section II-D, but specific details were not included. The 
purpose of this Appendix is to expand on the unusual concepts and 
geometry utilized in programming. It was not designed as a step by 
step analysis of the computer program nor as a discussion of concepts 
normally utilized in programming. However, for reference, a general 
list of arguments and arrays used is listed at the end of the Appendix. 

The two concepts to be discussed are the establishment of the jet 
oriented grid system and the computing of wind camponents parallel and 
nomal to the jet axis. In developing the concepts, for camputational 
purposes, the two and a half degrees latitude and longitude between 
grid points were normalized to one. Thus, the jet axis grid points 
could be taken off the grid as if they were located on a normal 
rectangular x-y grid. Distance between grid points in the y-direction 
was 150 nautical miles. In the x-direction, distance was adjusted for 
each latitude by the formula: DX=DY *COS (PHI), with PHI equal to the 
latitude. 

To determine the location of grid points 2.5 and 5.0 degrees lati- 
tude normal from the axis, the following scheme was used. First, the 
tangent at the individual jet axis grid point was determined. This 
was done by finding the distance in the x and y direction between the 
two jet axis grid points adjacent to the one used for calculation (see 
Fig. A-l). Using a computer library program ATAN 2 (A,B), the angle 


T was computed (Fig. AOQl). The sine and cosine of T, when added or 
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subtracted, as required, to the jet axis grid point, gives the location 
of the grid points 2.5 and 5.0 degrees latitude on elther side of the 
axis. Repeating this operation at each grid point along the jet axis 
established the grid described in Section II-D. 

To calculate the wind components normal and parallel to the jet 
axis, the angle T was used in conjunction with the u and v total 
balanced wind components at the curvilinear grid points. To obtain the 
parallel component, u was multiplied by cos T and v was multiplied by 
Sine T. To get the normal component, u was multiplied by sin T and v 
was multiplied by sinT (Fig. A-2). 

To assist any research in the future that might use this program, 
the following list of arguments and arrays is included. The complete 
program may be obtained from Dr. Russell L. Elsberry, Meteorology 
PeeEtmerit, Naval Postgraduate School, Monterey, California 93940. 

NP = number of grid points along the jet axis. 


NC 


number of grid points along the jet axis minus one. 
Used as a counter in DO-LOOPS 

NF, ND, NL, or NR = number of levels in the vertical. 

I = number of grid points along the Z or S-axis. 

K = number of grid points along the X or S-axis. 

J = number of grid points along the Y or M-axis. 

L = number of grid points along the Y or M-axis. 


ITT 


initial temperature field in degrees absolute. 


UU = initial zonal component of the total balanced wind in n/sec. 


VW = initial meridonal component of the total balanced wind in m/sec. 


WW = initial vertical motion field in mb/sec times 10°. 
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GP 


geopotential height fields in meters. 

PT = potential temperature. 
AC in front of the above five values means that it was the value on 
the curvilinear grid. 

PERVU = wind component normal to the jet axis. 

PARVU = wind component parallel to the jet axis. 
Any texm wi.th EDY before it designates an eddy value. All terms such 
as Al, A2, d, G, etc., refer to alphabetic designation given to various 


terms in equations discussed in the main part of this report. 
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Fig. A-i Schematic showing how angle T is calculated 
and used to define grid points 
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Fig. A-2 Schematic showing how angle T is used 
to acquire perpendicular and parallel 
wind components 
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